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ABSTRACT 

We present a picture of star formation around the Hll region Sh2-235 (S235) based 
upon data on the spatial distribution of young stellar clusters and the distribution 
and kinematics of molecular gas around S235. We observed 13 CO (1-0) and CS (2-1) 
emission toward S235 with the Onsala Space Observatory 20-m telescope and anal- 
ysed the star density distribution with archival data from the 2MASS survey. Dense 
molecular gas forms a shell-like structure at the south-eastern part of S235. The young 
clusters found with 2MASS data are embedded in this shell. The positional relation- 
ship of the clusters, the molecular shell and the H II region indicates that expansion of 
S235 is responsible for the formation of the clusters. The gas distribution in the S235 
molecular complex is clumpy, which hampers interpretation exclusively on the basis 
of the morphology of the star forming region. We use data on kinematics of molecular 
gas to support the hypothesis of induced star formation, and distinguish three basic 
types of molecular gas components. The first type is primordial undisturbed gas of 
the giant molecular cloud, the second type is gas entrained in motion by expansion 
of the Hll region (this is where the embedded clusters were formed), and the third 
type is a fast-moving gas, which might have been accelerated by winds from the newly 
formed clusters. The clumpy distribution of molecular gas and its kinematics around 
the Hll region implies that the picture of triggered star formation around S235 can 
be a mixture of at least two possibilities: the "collect-and-collapse" scenario and the 
compression of pre-existing dense clumps by the shock wave. 
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1 INTRODUCTION 

A mechanism for sequential star formation, in particular the 
formation of OB subgroups i n mo lecular clouds, was pre- 
sented by Hi megreen fe Ladal(|l977l ). Since that time, it has 
been ext ensively studied both theor etically and numerically 
(see e.g. lElmegreen &; Palousl (|2007h for recent results). The 
essence of the triggered star formation is the dynamical in- 
fluence of shock waves and winds from young massive stars 
and supernova explosions on their surroundings. The en- 
ergy input leads to shock front propagation and dense shell 
formation around massive stars or their clusters. A shock 
front can compress pre-existing dense clumps forcing them 
to contract and form new stars. The swept-up shell itself 
may become gravitationally unstable and fragment into new 
protostars. 

In their numerical study of star formation, triggered by 
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expanding Hll regions, iHosokawa &; Inutsukal (|20Q6h have 
argued that any expanding H n region is able to trigger star 
formation if the ambient molecular material is de nse enough. 
The s imilar conclusion has been reached by iDale et all 
(2007) on the basis of SPH modelling. However, these theo- 
retical statements require an observational check, especially 
in cases where the environment of an H n region possesses 
pronounced dumpiness and inhomogeneity. 

In this paper we present a study of the molecular envi- 
ronment of the H II region S235 (henceforth the S235 molecu- 
lar complex). The com plex is situated at a dis t ance between 
1.6 k pc and 2.5 kpc dGeorgelin et al.1 Il973l : llsrael fe Fell] 
1 19781 ). in the Perseus Spiral Arm. It appears to be a part 
of the giant molecular cloud G 174+ 2. 5, which has been ob- 
served in 12 CO (1-0) and 13 CO (1-0) lines bv iHever et all 
(1996). The S235 region has been created by the massive 
09.5V star BD +35° 1201 (|Georgelin et al.lll973T ). which has 
ionised and dispersed the surrounding molecular gas. The 
giant molecular cloud G174+2.5 harbors 3 other optical Hll 
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regions, namely, S231, S232, and S233. Shock waves, ex- 
panding away from these regions, should strongly influence 
surrounding molecular gas, changing significantly the kine- 
matic structure of the cloud. The S235 molecular complex 
is a site of active star formation, and objects at various 
stages of this process are observed there. First, there are 
young stellar clusters around S235. The clusters have been 
dete cted via analysis of the Spitzer data in vicinity of S235 
by lAllen et all (2005) and via inspection of 2MASS images 
toward IRAS source 05382+3547 bv lKumar et~ai1 (|2006l ). In 
the Section T3.2I of this paper on the basis of more extensive 
star counts we confirm presence of these clusters and pro- 
vide evidences for existence of some other infrared clusters 
in the region. Some members of these clusters are visible in 
POSS2/UKSTU Red and Infrared images taken from the 
Digitized Sky Survey archive. More deeply embedded and 
younger stellar/protostellar sources are present in the re- 
gion: the source IRAS 05382+354 7 was selected as a c andi- 
date precursor of ucHII region by Molina ri et al . (1991) on 
the basis of its colours. A 6.7 GHz methanol maser d etected 
toward this infrared source by Szymczak et al.l (2000) clearly 
indicates presence of an embedded young massive star. Note- 
worthy, IRAS 05382+3547 is situated furth er away from 
S235 H II region than the young stellar cluster ([Kumar et al.l 
2006). Our study shows that the separation between the cen- 
ter of the closest young cluster S235 Eastl and the IRAS 
source is considerable and amounts to about (Aa,AS)=(- 
80", 50"). Young stellar sources associated with methanol 
masers are usually f ound to be members o f stellar groups 
(jHunter et al.1 120061 : ILongmore et al.l l2006h . Consequently, 
the existence of a much younger stellar cluster, not yet pro- 
nounced even in Spitzer images, can be inferred. The pres- 
ence of various evolutionary stages suggests that star forma- 
tion in this region can be a sequential process. 



In this study we consider star formation scenario in the 
vicinity of S235. Specifically, we examine if the expansion of 
S235 is responsible for the birth of young clusters around 
it. We found specific features of star formation triggered by 
the expansion of an H n region, which were described in a se- 
ries of papers of Deharveng, Zay agno et al. (|Deharveng et al.l 
2003, 2005; Z avagno et all 2006): 1) radio emission in molec- 
ular lines as well as IR continuum emission of PAHs and dust 
show that there is a shell- like structure of molecular mate- 
rial close to the edge of the H n region and 2) positions of 
the young stellar clusters coincide with location of this gas 
shell. The procedure and results of star counts are described 
in Section 12.21 and Section 13.21 correspondingly. Observa- 
tions of emission in the CS (2-1) and 13 CO (1-0) lines are 
summarised in Section |2J"1 The spatial distribution and the 
kinematics of molecular gas are described in Sections 13.11 
and 13.31 We investigated kinematics of the dense molecu- 
lar gas (observed by CS (2-1) emission) and more rarefied 
gas (by 13 CO (1-0) emission) in details and contemplated 
that the observed kinematics might be one of the key instru- 
ments to distinguish between past, present and future stages 
of the triggered star formation process. This idea as well as 
the overall picture of star formation in the S235 molecular 
complex is presented in Section [4] 



Table 1. Basic observational parameters 
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13 CO(1-0) 


110201.353 


34 


0.068 


0.8 


CS (2-1) 


97980.968 


38 


0.076 


0.3 



2 OBSERVATIONAL DATA 

2.1 Molecular Emission 

We mapped the emission in the CS (2-1) and 13 CO (1-0) 
transitions toward the south-eastern part of S235 with the 
Onsala Space Observatory 20-m telescope in December 2005 
and February 2006. The observations in December 2005 
were dedicated to a search for peaks in emission in the 
CS (2-1) and 13 CO (1-0) lines around the 6.7 GHz methanol 
maser at a(J2000.0) = 05 h 41 m 33!80 s and £(J2000.0) = 
+35° 48 27" 00. We did extended mapping and observed 
emission in C 18 (1-0) and C 34 S (2-1) transitions toward 
selected positions in February 2006. We spent about 27 
hours of "ON" observational time for this project in Febru- 
ary 2006. The main observational parameters are presented 
in Table [1] The frequencies of observed em ission lines were 
taken from Lovas' catalogue (|Lovaslll992h . The full width 
at half maximum (FWHM) beam-size at both frequencies 
was about 40", and we used this value as a spacing for 
the CS (2-1) and 13 CO (1-0) maps. The main beam efficien- 
cies were about 0.54 for 97 GHz and 0.50 for 110 GHz. We 
used a bandwidth of 40 MHz with 1600 frequency channels, 
which corresponds to a velocity range of about 110 kms -1 
and a resolution of about 0.07 kms -1 . After the Decem- 
ber 2005 session we found that the emission in 13 CO (1-0) 
and CS (2-1) is extended, and performed observations in a 
frequency-switch mode. The system temperature, including 
the atmosphere, was typically 300-500 K for observations of 
CS(2-1) and 600-800 K for 13 CO(1-0). The typical noise 
level in maximum spectral resolution is 0.3 K for CS (2-1) 
and 0.8 K for 13 CO (1-0). The averaged uncertainties in the 
pointing accuracy were estimated to be 2.3" in azimuth and 
4.2" in elevation. 

2.2 Star Counts 

In order to find young star clusters we calculated the surface 
density distribution of the sources from 2MASS Point Source 
Catalog. This catalogue contains all stars detected in at least 
one of 3 filters (J, H and K s ) used for the survey. Every 
star which appeared in the catalogue was used for the star 
counts. Having made star counts at infrared wavelengths, 
we are able to find clusters which are not yet seen at visible 
wavelengths. Values of stellar density were deter mined at a 
recta ngular grid by the kernel estimator method ([Silver man! 
1986). We performed counts with different angular scales. 
The first map is a large-scale distribution of stellar density 
for a field of 100' x 100'. A grid of 201 x 201 nodes was used 
with a grid step of 0.5'. The kernel halfwidth was chosen to 
be l'. The second map was constructed for the vicinity of 
S235 with the dimensions of 20' x 20'. We used a grid of 
451 x 451 nodes, a grid step of 0.2', and a kernel halfwidth 
of 0.4'. Results are presented in Section [3721 
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3 RESULTS 

3.1 Spatial distribution of molecular gas 

Maps of CS (2-1) and 13 C0 (1-0) emission toward S235 are 
presented in Fig. [TJl. The origin of the map is a(J2000.0) = 
05 h 41 m 33!80 s and £(J2000.0) = +35°48 / 27 / /00. Multiple 
emission peaks are evident in both transitions. There is a 
peak in both transitions near (-40", 40"), while western 
peaks (Aa w —400") do not coincide. Regions with the most 
intense 13 CO (1-0) and CS (2-1) emission form a shell-like 
structure. Comparison of the maps wi th the radio c ontin - 
uum image of S235 at A21 cm made bv llsrael fe Fellil (|l978h 
shows that the shell surrounds S235 at the south-eastern 
side. This is also clearly seen in Fig. [TJ), which is a compos- 
ite of an optical image taken from the Digitized Sky Survey 
(POSS2/UKSTU Red) and the map of 13 CO (1-0) emission. 
At the inner side, this molecular shell is bordered by the dark 
filaments (white in the negative image) located to the south 
and south-east of S235 (Fig. [I})), which might well represent 
the dust collected by the shock wave. We suggest that the 
concentric dust and molecular shells have been shovelled by 
the expansion of the S235 Hn region. This suggestion will 
be further supported by consideration of the kinematics of 
the molecular and ionised gas. 

In Fig. [it we show the map of CS (2-1) emission su- 
perimposed on the negative 8.0/im image of the molecular 
complex taken from the Spitzer Space Telescope archive. 
The mid-IR emiss ion comes from th e strong PAH features 
at 7.7 and 8.6/xm ([Li &; Draind l2001). It is clearly seen that 
the dense clumps are bordered by the mid-IR patches from 
the side facing S235. This infrared emission appears as a re- 
sult of excitation of the PAHs in the thin layer at the border 
of the dusty region by the UV light of a close bright star or 
stellar cluster (e. g. IPovich et alJ (|2QQ7h ). 



3.2 Young stellar clusters in vicinity of S235 

We searched the environment of S235 for infrared star clus- 
ters to locate sites of recent star for mation in a much - 
wider area than in previous studies ([Allen et al ] 120051 : 
iKumar et al.ll2006h . The stellar density distribution in the 
region around S235 is shown in Fig. [2] Fig. [2^ shows 
the large-scale map (100' x 100'). Seven stellar clusters 
are found, listed with their coordinates and designations 
in Table [2] The richest cluster in the map is associated 
with the well- studied star- f orming region con taining objects 
S235A, S235B and S235C (jFelli et alJll997h . We designate 
it S235 A-B-C. Clusters S235 A- B-C, S235 E a stl an d Cen- 
tral are shown in the figures by lAllen et al. I (I2005T). Clus- 
ters S 235 Eastl and East2 were considered bv lKumar et al.l 
(2006). Two other clusters are situated more than 20' 
away from S235: these are the quite well studied bi- 
nary cluster S233IR between S231 and S233 (jPorras et al.l 
2000) and the relatively poor stellar condensation at the 
western edge of S232. There is also a cluster associated 
with IRAS source 5 361+ 3539 (G173.58+2.45) studied by 
IShepherd fe Watsonl <|2002h which is relatively poor and in- 
distinguishable among random condensations of the 2MASS 
sources. 

So, our analysis brought identification of 2 new clusters, 
S235 North- West and S232 IR (Fig. (St)). These clusters are 



Table 2. Coordinates of the young stellar clusters in the S235 
molecular complex 



Name 


R.A. (J2000.0) 


Dec. (J2000.0) 




(h m 


(° ' ") 


S235 East 1 


5 41 29.55 


35 49 21.65 


S235 East 2 


5 41 24.19 


35 52 15.96 


S235 Central 


5 41 07.85 


35 49 16.80 


S235 North- West 


5 40 45.73 


35 55 05.72 


S232 IR 


5 41 07.80 


36 12 27.00 


S233 IR 


5 39 10.95 


35 45 18.83 


S235 A-B-C 


5 40 51.94 


35 41 35.49 


IRAS 05361+3539 


5 40 54.45 


35 44 17.39 



relatively poor and we did not map them in CS (2-1) and 
13 CO (1-0) lines. 

Despite the thorough analysis, no other mid-IR clusters 
have been found in the 100' x 100' area around S235. It 
is noteworthy that five of seven bright clusters are located 
around S235 within an area with a radius of about 10'. So, 
S235 seems to be a focus for young clusters in the whole 
giant molecular cloud. 

The close vicinity of S235 was investigated with finer 
resolution. Contours of stellar density distribution, super- 
imposed on the map of CS (2-1) emission, are shown in 
Fig. EJ). It is clearly seen that clusters S235 Eastl, East2, 
Central, North- West and S235A-B-C consist of sub-clusters 
and there are small satellite groupings ar ound them. The 
exist ence of sub-clusters was reported by lAarseth &; Hillsl 
(1972) as a feature inherent to the young age of clusters, 
which is supposed to be less than one dynamical time. Our 
molecular line map covers the region containing three young 
clusters, S235 Eastl, East2 and Central. The locations of 
CS (2-1) clumps correlate with these clusters, which con- 
firms the extinction-based conclusion that they are still em- 
bedded in dense clumps of parental material. The densest 
parts of each cluster are projected onto the dense clumps. 
There is some displacement of CS (2-1) emission peaks rela- 
tive to peaks of the stellar density in the East2 and Central 
clusters which will be discussed in the next section. So, our 
data shows that young stellar clusters, associated with dense 
clumps of parental material, are embedded in the molecular 
shell-like structure around the H 11 region. 

3.3 Kinematics of molecular gas in the S235 
complex 

The locations of the young stellar clusters suggest that their 
formation might have been triggered by expansion of S235. 
As the vicinity of S235 is quite clumpy, conclusions based 
on morphology should be corroborated by the studies of 
gas kinemat i cs. Ki nematics of ionised gas was considered by 
lLafon et al.l (|l983h . He observed a velocity gradient from 
about -13 kms -1 to -30 kins" 1 from the South-East to the 
North- West. The data is incomplete and does not allow to 
study kinematics of the ionized gas in necessary details. So, 
systemic velocity of the molecular gas of the giant molecu- 
lar cloud G 174+ 2. 5 as a whole is more practical for fiducial 
point of the gas velocities. We discuss the velocity field of 
ionised gas in Sec. [4] in available details. Our observations 
show that in vicinity of S235 the majority of the 13 CO (1- 
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Figure 1. (a) Maps of 13 CO (1-0) (grey scale) and CS (2-1) (contours) emission toward S235. The contour levels for CS (2-1) map are 
shown from 1 K kms - 1 to 12 K kms -1 with a step 3 K kms -1 . The contour levels for 13 CO (1-0) map are shown from 5 K kms~ 1 to 
55 K kms -1 with a step 10 K kms -1 . Position of IRAS 05382+3547 and BD35°1201 are shown by asterisk and starlet correspondingly. 
The origin of the maps is a(J2000.0) = 05 h 41 m 33!80 s and 5(J2000.0) = +35°48 / 27 / /00. (b) Optical image of S235 taken from DSS 
(POSS2 Red filter) and overlayed map of 13 CO (1-0) emission (contours) from S235. Origin of the plot and the contour levels are the 
same as in panel (a). Positions of IRAS 05382+3547 and BD35°1201 are shown as in panel (a), (c) 8.0/im image of S325 from the 
GLIMPSE survey made with IRAC instrument from the Spitzer telescope archive and overlayed map of CS (2-1) emission (contours). 
Origin of the plot and the contour levels are the same as in panel (a). Positions of IRAS 05382+3547 and BD35°1201 are shown as in 
panel (a). 
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Figure 2. a) Maps of the stellar density around S235. The origin of the plot is the same as in Fig. [T] Contours correspond to a stellar 
density range from 13 stars • arcmin -2 to 31 stars • arcmin -2 with a step 2 stars • arcmin -2 . b) Map of CS(2-1) emission and stellar density 
distribution in the nearest vicinity of S235. The origin of the plot is the same as in Fig. ^ Contours correspond to a stellar density range 
from 20 stars • arcmin -2 to 56 stars • arcmin -2 with a step of 4 stars • arcmin -2 . Weak variations of background are omitted. Position of 
IRAS 05382+3547 is shown by asterisk. 
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0) and CS (2-1) emission has velocities in the range from 
-25 kms -1 to -15 kms -1 which is blue-shifted with respect 
to the overall velocity distribution of the giant molecular 
cloud which peaks at abou t -17,5 kms" 1 and has FWHP 
about 7 kms -1 (see Fig. 5 of iHever et al.l (|l996h ). Radial ve- 
locity of the ionising star BD +35° 1201 is about -18 kms -1 
as measured by the study of the photospheric absorption 
lines in the red optical range using 6-m BTA telescope at 
the Special Astrophysical Observatory of the RAS (Mlodik, 
Yushkin & Sobolev, private communication). 

Analysis of our observational data reveals three distinct 
kinematical components. Spatial distributions of 13 CO (1-0) 
and CS (2-1) emission in velocity ranges from -15 kms -1 to 
-18 kms -1 ("red" range), -18 kms -1 to -21 kms -1 ("cen- 
tral"), and -21 kms -1 to -25 kms -1 ("blue") are shown in 
Fig. [3^l and[3jl, ^jp and [3^, |3t and[3f, respectively. 

Gas responsible for emission in the "red" velocity range 
looks like three separate clumps visible in the 13 CO (1-0) 
line and situated at the edges of the map. There is no 
CS (2-1) emission (tracing dense gas) in the "red" range 
in the vicinity of S235, and there are no enhancements of 
stellar density toward these regions. In the large-scale map 
of 13 CO (1 -0) emission from the entire G174+2.5 molecu- 
lar cloud (|Hever et al.lll996h the emission, confined within 
the "red" velocity range, while also clumpy, appears to be 
distributed much more uniformly than the 13 CO (1-0) emis- 
sion in the other ranges of VW. Thus, we suggest that the 
gas emitting in the "red" range in our map represents qui- 
escent gas of the parental G 174+ 2. 5 molecular cloud, which 
was not disturbed by the expansion of S235 and can be con- 
sidered as "primordial" with respect to the gas where the 
young stellar clusters have been formed. 

Emission in the "central" velocity range within the re- 
gion, mapped in 13 CO (1-0) and CS(2-1), is spread more 
uniformly in our map in comparison with the other veloc- 
ity ranges. The 13 CO (1-0) emission with the "central" ve- 
locities is present almost in all spectra, except for the most 
south-western part of the map. This is shown in Fig. 2] which 
contains spectra of the 13 CO (1-0) line toward all observed 
positions. The transition from emission in predominantly 
"central" range to emission in predominantly "red" range is 
seen toward positions with double peaked 13 CO (1-0) lines. 

Emission from dense gas in the CS (2-1) line is observed 
mostly in the "central" velocity range, similarly to emis- 
sion in 13 CO (1-0) line (see Fig. and c). CS (2-1) emis- 
sion is clearly associated with the young stellar clusters (see 
Fig. [2J)). This is consistent with the assumption that the 
"central" velocity range represents gas, which has been af- 
fected and entrained in motion by the shock, propagating 
away from S235, and, thus, attained the velocity difference. 
Later, the young clusters were formed in the dense conden- 
sations of this gas. 

13 CO (1-0) and CS (2-1) emission in the "blue" range 
of Visr is confined to two well-defined patches which are lo- 
cated close to the S235 East2 and Central clusters. These 
are the most blue-shifted V\ sr values for the whole molecu- 
lar cloud G 174+ 2. 5, and this emission is observed only to- 
wards these two young clusters. The well-defined location of 
the "blue" clumps and the noticeable difference between the 
features seen in the "central" and "blue" ranges lead to the 
assumption that the "blue" patches might represent dense 



(CS (2-1)) and surrounding more rarefied ( 13 CO (1-0)) gas 
blown out by the stellar wind from the newly born clusters. 

The young clusters might have started to emerge from 
the parental material, pushing out surrounding molecular 
gas. The process wa s studied theoretically by, for example, 
IWeaver et all (j 19771 ). This process might produce a measur- 
able velocity difference between the dense gas moved by the 
stellar wind and the surrounding rarefied gas, which has not 
yet been impacted by the wind. Images of S235 East2 and 
Central clusters made in the K s -band (2MASS Survey) and 
at 4.5/xm (IRAC camera) contain filaments of dust which, 
according to their shapes, are most probably produced by 
the shock wave expanding from S235. Spectra of stars in the 
young clusters are needed to validate the idea about the ori- 
gin of the "blue" component. Indication that the stars have 
strong enough winds would be a reliable confirmation. 

Yet another kinematical evidence of the gas flows in 
the region is shown in Fig. [5] Profiles of 13 CO (1-0) and 
CS (2-1) lines toward the embedded clusters demonstrate 
relative velocity shift. The spectra toward ((Aa, AS) = (- 
40", 40")) have very good signal-to- noise ratio and we illus- 
trate the shift using this particular location (which is in the 
S235 Eastl cluster). We scaled the 13 CO (1-0) spectrum to 
the CS (2-1) line amplitude (see Fig. [5ji) and made a veloc- 
ity channel- by-channel subtraction of the CS (2-1) spectrum 
from it (see Fig. It is clear that the difference is large 
with respect to the noise level. A possible explanation of 
this velocity shift is the following: a wind blowing towards 
us compresses the gas and imparts on it a more negative V\ sr 
with respect to the undisturbed gas, which is uncompressed 
and, hence, relatively less dense. Since the CS (2-1) emis- 
sion traces more dense (compressed) gas than the 13 CO (1- 
0) emission (tracing undisturbed gas), we see the observed 
relative velocity shift. Note that both the 13 CO (1-0) and 
CS (2-1) lines peaks inside the same ("central") velocity 
range here. The non-detection of emission in the "blue" 
range toward the Eastl cluster might be related to a less 
evolved stage of the wind-driven outflow than in the East2 
and Central clusters. Probably, the wind from the Eastl 
cluster did not disperse the dense outflowing material yet, 
and we observe the velocity difference. 



4 DISCUSSION 

4.1 New view on star formation in the S235 
molecular complex 

The possibility of triggered star formation in the S235 molec- 
ular complex by expansion of the H n region has not been 
considered before. We found evidence supporting our view 
on star formation in this re gion in the Ha line emission data 
around S235 presented bv lLafon et al.l fl983). They found 
a gradient of Vi sr (Ha) from VW > -15 kms -1 in the south- 
ern part to V\sr < ~25 kms -1 in the northern part of S235 
where molecular gas is almost absent. Our map of 13 CO (1- 
0) emission does not contain data on this part of the S235 
complex, but channel maps of 12 CO (1 -0) and 13 CO (1-0 ) 
emission in this region can be found in IHever et al ] (|l996h . 
The channel maps allow suggesting that the northern part 
of the ionised nebula is expanding freely toward the ob- 
server. Data on radial velocities of radio recombination lines 
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Figure 3. Maps of 13 CO (1-0) (top row) and CS (2-1) (bottom row) integrated emission in different ranges of radial velocity superimposed 
on the optical image of S235 taken from the DSS. Origin of the plot is the same as in Fig.[T] (a) and (d) - 18 kms _1 < V\ sr < - 15 kms -1 , 
(b) and (e) -21 kms _1 < V lsr < -18 kms -1 , (c) and (f) -25 kms - ^ Vi sr < -21 kms" 1 . Contour levels are given for 10, 25, 50 and 70% 
from maximum values of the integrated intensity, which are 25 K kms -1 for 13 CO (1-0) and 10 K kms -1 for CS (2-1) emission. 



H109a bv lKazes. Walmslev fe Churchwelll ([19771 ) , on H167a 
and H140q bylsjjy erglate & Terzian (jl978l ) and on H91a by 



Qui reza et al. (2006)" toward some posit ions in S 235 are in 
agreement with the gradient found bv lLafon et al . (1983). 

The kinematic patt e rn of molecular gas was inter- 
preted bv lEvans fe Blair] (|l98lh (observations of 12 CO (1- 
0), 13 CO(1-0), C 18 O(l-0 ), H 2 CO(22 i -ln), HCN(l-O), 



HCO + (1-0) emission) and lLafon et all (|1983h (HCO + (l- 
0) and 13 CO (1-0) transitions) as a signature of existence 
of two molecular clouds with the mean radial velocities of 
-17 kms -1 and -20 kms -1 . Molecular gas to the south of 
S235, which contains the star forming regions S235A, S235B, 
and S235C, was interpreted as the "-17 kms -1 " cloud and 
gas to the east of S235, containing you ng clusters S235 East! 
and East2, as the "-20 kms -1 " cloud. lEvans fe Blairl (|l98lh 
mentioned that it is difficult to distinguish between the two 
clouds in directions where both velocity components are 
present in their map, but a thorough analysis of molecular 
kinematics was not done. 

Values of VW(Ha) around -17 kms -1 are observed in 
the direction of the southern part of the Hn region which 
is shielded by molecular gas around the S235 Central clus- 
ter. It is easy to see in Fig. [3] that the V\ ST of 13 CO (1- 
0) emission there is mainly in the range from -18 kms -1 
to -21 kms -1 , and absent at -17 kms -1 . It seems that 
molecular gas in front of the Hn region moves toward us 
faster than the ionised gas. This is expected if the gas is 
pushed b y the shock wave. It was sho wn b y numerica l mod - 
eling bv lHosokawa & Inutsukal (2006) and iHosokawal (2007) 
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Figure 5. Molecular spectra toward position with the offset coor- 
dinates (-40", 40"). a) CS(2-1) (grey) and scaled 13 CO (l-0)(l-0) 
(black) spectra, b) Residuals of spectra shown in the panel a). 



that the velocity difference between the shock and ionisa- 
tion fronts can be several kilometers per second if an Hll 
region is surrounded by uniformly distributed gas. The dif- 
ference becomes smaller if the gas density falls off with the 
distance from the H ll region as n ~ r~ w . It was shown by 
iFranco et al.l (|l990h that the critical value of w, when the 
initial ionised region still has a finite size, depends on a ratio 
of model parameters describing expansion of the H n regions 
in a non-uniform medium, w c ~ (1 — (r c /i?,s) 3 ) -1 , where r c 
is a radius of inner flat part in the density distribution and 
Rs is the Stromgren radius. We see that the molecular gas 
affected by the shock front (i.e., emitting in the "central" 
velocity range) propagates faster than the ionisation front, 
which means that the initial gas density distribution was 
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Figure 4. Spectral map of 13 CO (1-0) emission. Origin of the map is the same as in Fig. [T] Size of each cell equals 40" x40". The 
spectra are not smoothed. Scales of Vi sr and main beam temperature are shown in inset spectra ((Aq, A6)= (-80", 40")). 



not steep enough for formation of an unbound Hn region. 
Numerical modelling of the coupled chemistry and dynam- 
ics of an H II region, expanding into a clumpy medium with 
regions of decreasing and, probably, increasing gas density, 
is needed to clarify the history of star formation in S235 
molecular complex. 

The star formation in S235A-B-C is likely to be ini- 
tiated by a process different from the expansion of S235. 
These regions are situated quite far from the border of 
the S235, wh ich is delinea t ed in the radio continuum ob- 
ser yations by Israel & Felli (|l978l ) . Studies o f these regions 
bv lNakano fe Yoshidal (|l986h and lFelli et al . (1997) showed 
that the bulk of dense gas emission in the optically thin 
C 18 (1-0) and C 34 S (2-1) transitions appear in the VW 
range from about -18 kms -1 to -15 kins" 1 . This indicates 



that these star-forming regions are likely embedded in pri- 
mordial gas of the G 174+ 2. 5 giant molecular cloud. 



4.2 Scenario of star formation in the S235 
molecular complex 

There are several known ways for a n H II region to trig- 
ger star formation around itself (see lElmegreenl ([ 19981 ) for 
a review). In particular, contraction of pre-existing molec- 
ular clumps can be stimulated by the shock wave propa- 
gating away fro m the Hll region. This process was inves- 
tigated, e.g., by IBossI (|l995h for shocks with velocities of 
about 10 kms -1 an d low- mass molecular clouds, and by 
iKrebs &; Hillebrandtl |l983) for rapid shocks (faster than 
200 kms -1 ) and high -mass clouds. This scenario predicts 
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random distribution of dense clumps with embedded stellar 
clusters around an H II region. Ionised boundar y layers (sim- 
ilar to observed, e.g. bv lThompson et aT] [2004) are expected 
at the sides of the clumps which face the ionising star. 

Another way is the formati on of star-forming cores via 
"collect-and-collapse" scenario (jElmegreen fe LadJl977h . It 
predicts the existence of young clusters embedded in the 
shell of molecular gas located ahead of the ionisation front. 
In an ideal spherically-symmetric case this gas is expected 
to be observed as a ring of molecular emission around an 
Hn region. There should also be rings of mid-IR emission 
produced by PAHs in the PDR region and millimetre dust 
emission around the Hn region. The young clusters have 
to be situated on a neutral side of the compressed layer, 
between the PDR and the molecular gas of the parental 
cloud. 

There are three morphological signatures of the "collect- 
and-collapse" scenario here: 1) young stellar clusters are em- 
bedded in the dense molecular shell around S235, 2) PAHs 
residing in the PDR (see Fig. [I]: for emission at 8.0/xm) pro- 
duce the shell-like mid-IR emission, 3) the young clusters 
S235 Eastl and East 2 are situated on the neutral side of 
the mid-IR shell. The mid-IR shell looks broken toward the 
S235 Central cluster, which is projected on the optical im- 
age of the Hll region, but Eastl and East2 are located to 
the east of S235. All three clusters are embedded in clumps 
emitting in the "central" velocity range. This fact is an evi- 
dence supporting ongoing "collect-and-collapse" scenario. 

It is not easy to single out one scenario of triggering 
of star formation by S235 because the molecular environ- 
ment of the Hll region is clumpy. The actual scenario can 
be a mixture of the scenarios mentioned before. The answer 
can be prompted by consideration of evolutionary stages of 
the embedded clusters. The age similarity between the three 
clusters would rather be an argument in favour of a "collect- 
and-collapse" scenario, as it seems less probable to find three 
independent pre-existing clumps with similar mass and den- 
sity on a similar distance from the H n region. On the other 
hand, even in a "collect-and-collapse" scenario, clusters do 
not need to be exactly coeval, the age distribution might be 
caused by different densities of initial clumps in a fragment- 
ing molecular shell. We discuss evolutionary stages of the 
embedded clusters below. 

4.3 Evolution stages of the young embedded 
stellar clusters 

Our data on the distribution of the stellar density and the 
dense gas as well as the kinematic data indicate that the 
clusters S235 East2 and Central are probably more evolved 
than Eastl. Peaks of the stellar density are displaced relative 
to the peaks of CS (2-1) emission for the S235 Central and 
East2 clusters. These clusters might have already started ac- 
tive dispersion of the dense gas by their stellar winds, giving 
rise to the "blue" component of the 13 CO (1-0) and CS (2- 
1) clumps in Fi g . E3s. Obser vations of Herbig Ae/Be stars 
bv iFuente et~aD (|l998l 2002) showed a sequence of objects 
which represents a history of mass dispersal around them: 
from embedded Herbig Ae/Be stars to stars located in a 
cavity within the molecular cloud. They suggested that en- 
ergetic winds from the stars swept out gas and dust, forming 
a cavity. Similar phenomenon can take place in S235. 



In contrast, the S235 Eastl cluster does not show notice- 
able displacement of the molecular emission peak relatively 
to the stellar density maximum and does not show any asso- 
ciated "blue" range emission. This cluster might be younger. 
Alternatively, S235 Eastl can be at the same evolutionary 
stage as the other two clusters, if it does not contain stars 
which are luminous enough to produce a fast and powerful 
impact on its environment. Data on C 18 (1-0) emission to- 
ward the clusters are needed to estimate masses of the three 
dense clumps and those of the gas emitting in the "blue" 
velocity range. 

The question about relative ages of the clusters could 
be answered by sp ectral studies of their stellar population. 
lAllen et a l. (2005) identified some of young stellar objects, 
belonging to the S235 Eastl, East2 and Central clusters, as 
Class I and Class II sources. They found Class I objects in 
all three clusters, but Class II objects were discovered only 
in the S235 East2 and Central clusters. This result is in 
agreement with our guesses described above. But this result 
cannot be considered as a strong evidence of age spread be- 
cause of small number of identified objects (two young stars 
in S235 Eastl, seven in S235 East2 and 10 in S235 Central). 
More objects must be measured to make definite conclu- 
sions. 

It is tempting to suggest that kinematic components, 
which we see in S235, correspond to sequential evolutionary 
stages. In this line of reasoning, the most advanced stage 
of the triggered star formation would be represented by the 
East 2 and Central clusters, where a non- negligible amount 
of gas has already been accelerated by stellar winds. The 
clump coincident with the Eastl cluster would be at the 
intermediate stage as it is actively forming stars but have 
not yet developed any observable wind signatures. 

Interesting is also the clumpy "red" component, consti- 
tuting the structure of the parent molecular cloud. We do 
not know for sure if the three "clumps" seen on the 13 CO (1- 
0) map are not just projection effects, but real physical enti- 
ties. But if they are real, they may represent a visible part of 
a more abundant population of clumps, waiting for some ex- 
ternal influence to trigger the formation of star clusters and 
massive stars in them. As such, these clumps are of particu- 
lar interest, and we plan to continue the observational study 
of the S235 complex, with the ultimate goal to check if we do 
see objects in this area which represent three distinct stages 
related to past (S235 East2 and Central), present (Eastl), 
and future (hypothetical "red" clumps) of the triggered for- 
mation of star clusters and massive stars. 



5 CONCLUSIONS 

We investigated the kinematics and distribution of molecular 
gas and the spatial distribution of young stellar clusters in 
the S235 molecular complex and found evidences that the 
process of star formation in the complex was triggered by 
expansion of the H n region. 

The H II region S235 is surrounded by a shell-like enve- 
lope with dense molecular clumps at the eastern and south- 
ern sides. The shell circles the ionisation front and is sep- 
arated from the Hll region by a PDR envelope visible in 
mid-IR. The dense molecular clumps contain young stellar 
clusters not yet seen at visible wavelengths. We suggest that 
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the expansion of S235 is responsible for the formation of the 
clusters. 

Emission from molecular gas in the 13 CO (1-0) and 
CS (2-1) transitions reveals three kinematically distinct 
components, which can be interpreted in terms of the 
triggered star formation scenario, presumably representing 
three stages in the influence of young stellar objects on the 
surrounding gas. The clumpy gas emitting in the "red" ve- 
locity range (-18km s _1 < Vi sr <-15kms _1 ) may represent 
primordial material of the G 174+ 2. 5 giant molecular cloud 
undisturbed by expansion of S235. This gas is not dense 
enough to be readily observable in CS (2-1) emission. Emis- 
sion in the "central" range (-21 krns" 1 < VW <-18 krns" 1 ) 
can be produced by the gas clumps which were affected 
and entrained in motion by the expansion of S235. These 
are the clumps where the young clusters were formed. Re- 
moval of molecular gas by the stellar winds from the clusters 
S235 East 2 and Central might have given rise to the "blue" 
component (-25 krns -1 < VW <-21 krns -1 ) in 13 CO (1-0) 
and CS (2-1) emission. 

The cluster S235 Eastl is probably younger and rep- 
resents an earlier evolutionary stage than the Central and 
East2 clusters. Another possibility is that S235 Eastl con- 
tains less luminous stars, being at a similar evolutionary 
stage. More extensive spectral data on the stellar popula- 
tion of these clusters are needed to determine their relative 
ages. 
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